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1.  Introduction 


Free-space  lasers  are  becoming  more  widely  used  in  the  military  environment,  and  it  is  important 
to  develop  a  capability  to  detect  and  remotely  characterize  laser  beams  from  an  off-axis  location. 
Many  applications  in  communication,  surveillance  and  targeting,  require  relatively  high-power 
or  pulsed  laser  sources  and  involve  propagation  distances  of  several  kilometers.  If  such  lasers  are 
used  as  a  threat,  it  would  be  desirable  to  know  the  power  and  location  of  the  source  and  the 
direction  of  the  beam  in  order  to  provide  warning  and  enable  counter-measures.  Previous  work 
has  demonstrated  that  pulsed  laser  beams  can  be  detected  from  scattering  by  aerosols  in  the  path 
and  the  beam  can  be  imaged  with  suitable  detectors  [1].  There  a  10x10  detector  array  was  used 
and  the  scattered  intensity  was  found  to  be  roughly  consistent  with  the  prediction  based  on  Mie 
scattering  from  a  generic  continental  aerosol  distribution.  Aerosol  concentrations  can  vary 
widely  depending  on  regional  and  local  weather  conditions.  More  recent  work  has  investigated 
how  details  of  the  aerosol  distribution  would  affect  the  scattered  intensity  [2,3].  For  example, 
increased  humidity  will  lead  to  larger  aerosols  and  generally  enhance  the  scattering.  In  typical 
conditions,  even  at  low  altitude,  the  scattering  is  extremely  weak  and  can  be  difficult  to  detect  in 
the  presence  of  background  light. 

In  this  work,  we  present  a  detailed  theory  of  the  remote  detection  of  scattered  light  and 
make  comparisons  with  Monte  Carlo  calculations  to  understand  the  effect  of  multiple  scattering 
on  the  received  intensity.  Extensive  experimental  intensity  and  temporal  data  from  a  high-speed 
30-channel  angle-resolved  receiver  is  presented  and  compared  with  the  analytic  theory.  Volume 
scattering  is  estimated  by  using  a  range  of  modeled  aerosol  distributions  representative  of  the 
observed  weather  conditions.  The  high  angular  resolution  of  the  receiver  allowed  the  direction  of 
the  source  to  be  tracked  as  the  geometry  of  the  experiment  evolved  over  time.  Assuming  a 
volume  scattering  function,  the  laser  power  and  distance  as  well  as  the  direction  of  the  beam  all 
affect  the  scattered  intensity  at  the  receiver  and  these  parameters  cannot  be  readily  distinguished 
from  intensity  data  alone.  If  the  location  of  the  source  is  know,  there  would  still  be  an  ambiguity 
between  power  and  beam  direction.  In  this  case,  we  show  that  accurate  timing  information  of  the 
received  scattered  light  can  be  used  to  estimate  beam  direction  and  subsequently  laser  power. 


2.  Off-Axis  Scattering  Theory 

The  scattering  plane  geometry  is  defined  in  Fig.  1 .  The  receiver  and  laser  source  are  separated  by 
R  and  the  beam  from  the  source  is  oriented  at  an  angle  t//  from  the  baseline  between  them.  At  the 
size-scale  of  interest  in  this  work  it  is  a  good  approximation  to  neglect  beam  size  and  divergence 
and  treat  the  beam  as  one-dimensional.  The  receiver  collects  scattered  light  from  a  laser  pulse  at 
an  angle  6  from  the  baseline  corresponding  to  a  propagation  distance  z  along  the  beam  line.  The 
scattered  light  travels  a  distance  r  to  the  receiver.  The  angle  from  the  direction  of  the  source  to  a 
reference  direction  on  the  receiver  is  given  by  &s. 

We  consider  the  light  collected  by  the  receiver  to  be  only  due  to  volume  scattering  from 
aerosols  in  the  path  of  the  beam.  The  scattering  is  characterized  by  a  volume  scattering  function 
fi(y/+9)  (m'sT1)  which  depends  on  the  scattering  angle  ^+6* from  the  direction  of  propagation.  In 
general,  the  scattering  will  vary  with  position  but  for  simplicity,  we  assume  /?  to  be  spatially 
uniform.  The  laser  power  is  extinguished  along  the  beam  line  as  Poexp(-az)  where  a  is  the  beam 
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extinction  coefficient  due  to  absorption  a  and  total  scattering  b  =  2/rJ()  sin (d)(%d)dd  along  the 

path.  We  will  neglect  multiple  scattering,  which  is  a  good  approximation  if  the  mean  free  path 
1  lb  is  large  compared  with  the  spatial  scale  of  the  experiments.  If  we  assume  a  uniform 
homogeneous  medium,  the  angle  and  time-dependent  scattered  radiance  at  the  receiver  from  a 
laser  pulse  with  an  initial  power  Ppt)  can  be  written  in  terms  of  the  beam  extinction  coefficient  a 
and  the  volume  scattering  function  f. 3  as 


dl(d,t) 

dO 


=  po  \_tsit  (0  - 1  ]exp[-a(z  +  r)] 


fXd+y/)  dz 

7  ~dd' 


(1) 


Here  z(0)=Rsm(d)/sm(d+y/)  and  r(6)  =  /?sin(t//)/sin(t9+  <//)  are  the  distances  from  the 
scattering  volume  to  the  source  and  receiver  respectively  and  (sr(  d)  =  (z+r)/c  is  the  propagation 
time  along  the  path  from  the  source  to  receiver.  If  a  detector  collects  scattered  light  over  a 
limited  one-dimensional  angular  range  Ad  in  the  scattering  plane,  the  total  intensity  /( d.t) 
centered  at  0  in  the  middle  of  the  range  A  6  will  be  a  convolution  of  the  time-dependent  pulse 
with  the  scattering  along  the  path  subtended  by  Ad.  After  simplifying,  this  becomes 


I(d,t)=  J  P0  [tSR id’)  -  /]exp  \—a(z  +  r)]  ^ .  -  ¥)  dd' .  (2) 

e-Ken  R  smW 


The  peak  intensity  at  the  receiver  is  usually  required  and  therefore  it  is  necessary  to 
determine  the  time  t  that  maximizes  Kd.t).  The  propagation  time  increases  with  both  y/  and  d  and 
the  variation  of  /sr  with  d  is 


" jp-  =  ^sin(^)/ cos[(<9+  y/)/2]2 . 


(3) 


Some  approximations  that  simplify  Eq.  2  can  often  be  made,  depending  on  the  variation  of 
P<)(tsR(d)-t)  and  /3{0±y/)  over  Ad  [1  J.  Consider  the  ratio  r/AtSR  of  laser  pulse  width  rto  the  range 
of  ?sr  over  Ad.  In  the  long-pulse  limit  where  the  ratio  is  large,  the  maximum  intensity  occurs  at  t 
~  tsR(d)  and  Po(tsR( d)-t)  ~  Pq(Q).  As  the  ratio  decreases  and  t  becomes  comparable  to  or  shorter 
than  AtSR,  the  sensitivity  of  Kd.t)  with  respect  to  t  increases.  The  maximum  intensity  is  found 
when  the  pulse  is  centered  at  angles  nearer  to  the  source  because  the  scattering  function  ft  is 
generally  larger  for  smaller  scattering  angles.  In  the  analysis  below,  the  integration  (Eq.  2)  was 
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replaced  by  an  analytic  summation  over  discrete  angles  within  Aft  and  the  sum  was  maximized 
with  respect  to  t  as  a  parameter.  A  much  simpler  and  often  good  approximation  is  to  simply 
center  the  pulse  at  the  midpoint  of  integration.  For  our  experimental  conditions  (see  below),  we 
found  this  approximation  to  underestimate  the  maximum  of  1(6, t )  by  less  than  -10%  for  larger  y/ 
and  smaller  ft  where  the  error  is  greatest. 


3.  Experiment 

Scattering  data  at  ft  =  1.06  pm  was  collected  with  the  Ship  Board  Laser  Acquisition  System 
(SBLAS)  [4],  a  high-speed  multi-channel  receiver  that  could  record  arrival  time  and  intensity 
over  a  horizontal  field  of  view  (FOV)  from  -45°  to  +45°.  The  SBLAS  was  composed  of  two  15- 
element  linear  detector  arrays  with  separate  collection  optics  that  covered  the  range,  in  receiver 
coordinates  6- 6s,  from  -45°  to  0°  and  from  0  to  45°  respectively.  The  aperture  normal  and  optical 
axis  of  each  side  was  inclined  to  the  middle  of  the  respective  45°-range.  Each  of  the  30  angular 
channels  had  a  narrow  asymmetric  FOV,  Aft  =  3°  in  the  horizontal  plane  and  -15°  <  6fy  <  15°  in 
the  vertical  plane.  Each  channel  digitized  the  signal  and  recorded  the  crossing  time  over  9 
threshold  values.  Depending  on  the  highest  threshold  crossed,  the  true  peak  intensity  could 
exceed  the  threshold  value  by  a  factor  of  one  to  four.  The  experiments  described  below  were 
conducted  in  full  daylight  conditions  in  the  presence  of  significant  background  radiance. 
However,  the  detectors  were  AC-coupled  and  the  constant  background  had  little  effect.  False 
alarms  in  the  receiver  were  further  reduced  by  requiring  a  valid  detection  event  to  consist  of 
responses  in  multiple  channels  with  timing  consistent  with  the  propagation  of  a  laser  pulse. 

Beam  scattering  experiments  were  performed  over  several  days  off  the  coast  of  Cowley 
beach,  Queensland  Australia.  A  laser  source  (ft  =  1.06  pm)  with  pulse  width  x  =  20  ns  full  width 
at  half  maximum  was  located  on  the  beach  about  1-2  m  above  sea  level  and  the  beam  was 
collimated  and  directed  at  a  stationary  target  boat  offshore.  The  SBLAS  receiver  was  mounted 
on  another  ship  with  its  axis  facing  to  starboard  and  perpendicular  to  the  heading  of  the  ship.  The 
extensive  experimental  data  taken  on  1  July  2004  exemplifies  the  capability  of  the  receiver  over 
a  significant  variation  of  scattering  geometry  and  allows  a  meaningful  comparison  with  the 
theoretical  model  discussed  above.  GPS  positioning  data  taken  on  the  ship  allowed  the  path  of 
the  receiver  relative  to  the  laser  beam  to  be  reconstructed  as  shown  in  Fig.  2.  During  this 
experiment  the  ship  followed  a  path  -  4.8  km  from  shore  and  approximately  perpendicular  to  the 
laser  beam. 

As  mentioned  earlier,  atmospheric  aerosols  in  the  beam  path  are  the  source  of  off-axis 
laser  scatter.  Some  limited  measurements  of  the  aerosol  size  distribution  at  the  ship  were  made 
however  they  were  not  considered  to  be  of  sufficient  quality  to  predict  the  volume  scattering 
accurately.  Furthermore,  over  the  beam  path  that  extended  nearly  5  km  from  the  shore,  the 
aerosol  distribution  was  unlikely  to  be  uniform.  Fortunately  there  has  been  significant  effort 
devoted  to  developing  models  that  describe  the  generation  and  distribution  of  aerosols  in 
maritime  environments  based  on  environmental  conditions  [5,6,7],  Here  we  have  used  the 
Advanced  Navy  Aerosol  Model  (AN AM)  [7]  to  predict  representative  aerosol  distributions 
consistent  with  the  weather  conditions  preceding  and  during  the  experiment.  ANAM  includes 
five  distinct  types  of  uniform  spherical  particles,  each  with  a  log-normal  distribution  based  on 
radius.  Details  of  the  model  are  described  in  the  appendix.  Over  open  water,  the  current  wind 
speed  and  the  preceding  24-hour  average  wind  speed  are  perhaps  the  most  significant  factors  in 
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the  production  of  aqueous  aerosols,  especially  the  larger  aerosols  that  lead  disproportionally  to 
small-angle  forward  scattering. 

Measurements  of  wind  speed  w,  temperature  T  and  relative  humidity  RH  were 
continuously  recorded  on  the  ship  beginning  at  9:00  AM  local  time.  The  collection  of  scattering 
data  began  at  approximately  2  PM  and  lasted  for  -280  seconds  while  the  ship  followed  the  path 
shown  in  Fig.  2.  During  this  period  the  wind  speed  varied  between  2  and  8  m/s  and  the  average 
wind  speed  over  5  hours  preceding  the  test  was  -3  m/s.  It  is  important  to  note  that  these 
conditions  were  measured  some  kilometers  from  the  beam  scattering  and  were  probably  not 
uniform  over  the  entire  region.  Therefore  we  used  ANAM  to  generate  four  representative  aerosol 
distributions  that  should  have  bracketed  the  range  of  scattering  conditions  during  the  tests.  The 
environmental  parameters  used  in  the  four  cases  are  given  in  Table  1.  The  largest  and  heaviest 
aerosols  with  modal  radius  -10  pm  are  generated  by  wind  action  on  the  sea  surface  and  the 
equilibrium  density  decreases  quickly  with  height  h.  They  are  expected  to  be  significant  near  the 
sea  surface  at  h  =  2  m.  For  comparison,  a  case  with  h  =  300  m  was  chosen  to  effectively 
eliminate  these  aerosols  which  were  not  present  in  the  earlier  NAM  aerosol  model  [5].  The  air- 
mass  parameter  (amp)  is  meant  to  describe  the  relative  coastal  influence  (from  0  to  10)  and 
lacking  specific  information,  we  chose  cases  with  amp  =  3  (less  coastal  influence)  and  amp  =  8 
(more  influence). 

The  ANAM  model  includes  wavelength-dependent  absorption  and  refractive  indices  for 
the  spherical  aerosols.  For  each  of  the  four  cases,  the  Bohren-Huffman  MIE  scattering  code  [8] 
was  used  to  calculate  the  angle-dependent  volume  scattering  function  appropriate  for  detection  of 
scattered  light  from  a  source  polarized  perpendicular  to  the  scattering  plane,  (5  =  (A,I2ti)2(S\\  - 
Su),  where  S,,  are  Mueller  matrix  elements.  These  representative  scattering  functions  for  A  =  1.06 
pm  are  compared  in  Fig.  3.  Some  notable  features  are  apparent.  First,  there  is  a  pronounced 
decrease  in  small-angle  forward  scattering  for  the  aerosols  of  Case  3  that  do  not  include  the 
largest  radius  aerosols  centered  near  10  pm.  Case  4  differs  from  Case  2  only  by  a  larger  air-mass 
parameter  and  consequently  the  inclusion  of  the  small,  high  refractive  index,  dust-like  aerosols 
that  lead  to  an  increase  in  scattering  above  about  3°.  The  total  beam  extinction  (a+b)  varied  over 
a  factor  of  -5  for  these  cases.  Measurements  of  the  laser  beam  extinction  would  have  been 
useful  to  help  characterize  the  aerosol  conditions  during  the  experiments  but  were  not  taken. 

Next  we  compare  the  model  predictions  with  the  data  from  the  SBLAS  receiver  during 
the  experiment  off  Cowley  Beach,  Australia.  GPS  position  data  taken  on  the  ship  were  recorded 
along  with  the  receiver  data  at  a  rate  of  10  readings  per  second  and  were  precise  to  about  2  m  and 
5  m  along  the  east-west  (x)  and  south-north  (y)  directions  respectively.  A  piece-wise  linear  fit  of 
the  GPS  coordinates,  in  segments  of  15  seconds,  was  used  to  construct  the  ship  heading  from 
which  the  source  angle  6s  or  orientation  of  the  receiver  relative  to  the  source-receiver  baseline 
and  beam  angle  y/  were  calculated.  In  Fig.  4  we  show  how  6s  fluctuated  2-3°  during  the 
experiment  due  to  variations  of  the  heading  from  a  straight  path.  However  the  beam  angle  y/, 
because  it  is  defined  relative  to  the  position  (rather  than  orientation)  of  the  receiver,  was 
insensitive  to  the  small  heading  fluctuations.  An  estimate  of  the  source  angle  6se  was  obtained 
from  the  intensity  data  by  taking  the  angle  of  the  channel  with  the  largest  signal.  Since  the 
direction  of  the  beam  to  the  left  or  right  was  always  apparent,  ambiguities  among  multiple 
channels  recording  the  same  intensity  threshold  were  resolved  by  selecting  the  channel  nearest 
the  source.  Although  6se  was  only  resolved  to  3°,  these  measurements  are  clearly  correlated  with 
the  GPS  tracking  data. 
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Scattering  data  from  the  SBLAS  receiver  was  recorded  for  each  firing  of  the  laser  and 
2742  measurements  were  collected  as  the  receiver  traversed  the  beam  line.  The  SBLAS  signal 
was  corrected  to  give  intensity  (W/m  )  in  the  direction  of  receiver  coordinate  0-  Os  and 
normalized  to  the  peak  power  of  the  laser  pulse,  assuming  a  Gaussian  temporal  profile  with  pulse 
width  t  =  20  ns  full  width  at  half  maximum.  Six  measurements  from  discrete  laser  firings  at 
selected  beam  angles  y/  in  the  range  from  -6.5°  to  +6.5°  are  shown  in  Figs.  5  and  6.  The  data 
generally  lie  within  the  range  of  predictions  of  the  scattering  model  (Eq.  2)  using  the  volume 
scattering  functions  associated  with  the  four  aerosol  cases.  In  particular,  the  data  show 
essentially  the  same  exponential  decrease  of  intensity  with  angle  over  most  of  the  angle  range.  At 
small  scattering  angles  0 4-  y/,  less  than  about  2°,  the  model  predicts  a  large  enhancement  because 
of  the  increase  of  both  (3(0+ y/)  and  dz/d#(Eq.  1).  This  is  especially  pronounced  when  the  largest 
aerosols  are  present  and  contribute  to  the  small-angle  volume  scattering  in  cases  1,2  and  4.  Some 
evidence  of  a  slight  enhancement  in  the  data  was  seen  for  I  yA  <  ~1°  but  generally  the 
measurements  do  not  show  the  same  magnitude  of  small-angle  scattering,  perhaps  in  part 
because  fewer  large  aerosols  were  actually  present  in  the  beam  path. 

Although  the  single-scattering  model  gives  quite  good  agreement  with  the  measured 
intensity,  it  is  important  to  understand  how  much  the  accuracy  is  affected  by  the  neglect  of 
multiple  scattering.  Since  each  angular  channel  within  the  SBLAS  receiver  has  a  large 
asymmetric  FOV,  some  multiple  scattering  especially  from  laser  power  scattered  close  to  the 
beam  line,  will  be  detected  and  lead  to  a  larger  signal  than  predicted.  In  order  to  verify  the  model 
we  have  performed  a  Monte  Carlo  (MC)  simulation  for  a  geometry  similar  to  the  experiment 
using  the  volume  scattering  function  of  Case  4  which  has  the  largest  beam  extinction 
(~  16%/km).  The  MC  procedure  has  been  described  in  previous  work  [9].  In  this  simulation  108 
photons  were  launched  along  a  beam  line  normal  to  a  collection  plane  located  d  =  4778  meters 
away.  The  normalized  scattered  intensity  (m~2)  for  a  desired  beam  angle  (//was  obtained  from  the 
fraction  collected  in  an  annular  ring  with  center  radius  x  =  d  tan(  y/)  and  width  Ax  =  5  m  and  then 
resolved  into  angular  channels  A0  =  3°.  A  scattering  plane  for  each  photon  was  defined  by  the 
beam  line  and  the  impact  point  in  the  annulus  which  then  allowed  the  out-of-plane  angle  of 
incidence  to  be  calculated.  Although  the  SBLAS  receiver  had  a  symmetric  vertical  FOV,  the 
proximity  of  the  ocean  surface  to  the  beam  and  receiver  meant  that  the  effective  vertical  FOV 
was  reduced  to  0  5  <15°  and  therefore  only  photons  collected  within  this  FOV  were  counted. 

Figure  7  compares  normalized  receiver  intensities  (W/m2/W)  for  beam  angles  y/  =  0.6°,  2.4°  and 
9.5°  (  0s  =-ys)  with  the  model  prediction  where  the  laser  power  P0  was  constant.  Both  the  model 
and  the  MC  results  show  the  enhanced  receiver  intensity  at  small  scattering  angles  0 4-  y/  and  the 
nearly  exponentially  decrease  at  larger  angles.  As  expected,  the  single-scatter  model  somewhat 
underestimates  the  intensity  when  the  beam  is  in  the  receiver  FOV.  The  MC  results  also  show 
weak  intensity  at  angles  behind  the  beam  due  to  multiple  scattering  that  is  not  accounted  for  in 
the  model.  When  only  singly-scattered  photons  were  counted,  the  MC  results  were  consistent 
with  the  model.  When  the  beam  was  in  the  receiver  FOV,  the  relative  difference  increases  with 
larger  beam  angles  because  each  channel  FOV  encompasses  a  larger  volume  of  space  from 
which  multiply  scattered  light  originates.  At  y/=  9.5°  the  difference  is  -15%,  however,  including 
scattering  from  a  symmetric  FOV  -15°  <  0y  <15°  would  have  doubled  this. 

One  would  like  to  characterize  the  power  and  location  of  the  laser  source  as  well  as  the 
direction  of  the  beam  from  the  receiver  data  and  the  model  above.  For  many  conditions, 
including  those  of  this  experiment,  the  analytic  model  provides  a  good  approximation  to  the 
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scattered  intensity  observed  remotely.  However,  assuming  that  the  volume  scattering  can  be 
estimated  accurately,  these  three  parameters  cannot  be  determined  unambiguously  from  the 
receiver  intensity  alone  (Eq.  2).  Two  or  more  spatially  separated  receivers  could  be  used  to 
locate  the  direction  and  distance  to  the  laser  source  by  triangulation  and,  in  this  case,  the  relative 
intensity  could  then  be  used  to  resolve  the  ambiguity  between  laser  power  and  beam  direction. 
Here  we  demonstrate  that,  with  a  single  receiver,  the  beam  angle  if/  can  be  determined  from 
timing  measurements  of  the  scattered  light  if  the  position  of  the  source  is  known.  Then,  with  this 
information,  the  laser  power  can  be  estimated. 

The  SBLAS  receiver  recorded  the  relative  times  if  and  when  each  of  the  level  thresholds 
was  crossed  for  each  receiver  angle.  We  considered  the  crossing  of  the  highest  threshold  fa  to 
approximate  the  time  of  peak  intensity  and  representative  timing  data  fa  over  the  receiver  angles 
with  some  detected  signal  for  selected  beam  angles  if/  are  shown  in  Fig.  8.  The  analytic  path  time 
tsR  ~  to  =  (z+r)/c  -  to,  based  on  the  known  geometry,  matches  this  data  with  an  offset  time  to 
obtained  by  a  least  squares  linear  fitting  procedure.  Here,  the  variation  of  /sr<  with  angle  is 
essentially  linear.  There  is  less  than  1%  variation  between  d/sk/d#  calculated  at  <0>  and 
<dfsR/d#>  averaged  over  6.  For  each  laser  shot,  the  experimental  average  time  derivative 
<dtn/d#>  was  calculated  from  a  linear  fit  of  fa  over  angle.  These  measured  values  agreed  well 
with  the  predicted  derivative  (Eq.  3),  with  devaluated  at  the  mid-point. 

The  beam  angle  if/c<m  be  determined  from  d/sk/dd  using  Eq.  3  if  the  source  angle  ds  and 
range  R  are  known.  For  each  of  the  2742  measurements,  an  estimate  of  if/  was  calculated  by 
substituting  <dfaldO>  on  the  left  hand  side  of  Eq.  3  and  solving  for  fusing  the  estimated  source 
angle  Ose  and  known  range  R  (from  GPS  data).  Figure  9  shows  the  difference  between  the 
calculated  beam  angle  y/c  and  the  angle  from  the  GPS  tracking  data.  Over  the  range  from  -6°  to 
+6°  the  rms  error  is  0.22°.  The  laser  power  can  then  be  calculated  from  Eq.  2  based  on  an 
estimate  of  the  beam  extinction. 

4.  Conclusion 

The  remote  detection  of  pulsed  laser  beams  by  scattering  from  atmospheric  aerosols  has  been 
investigated.  A  theory  of  the  angle  and  time-dependent  scattered  intensity  for  a  general 
configuration  of  laser  source,  receiver  and  beam  direction  has  been  described.  We  have 
presented  extensive  measurements  of  off-axis  scattering  observed  at  ranges  up  to  ~5  km  with  a 
high-speed  multi-channel  receiver.  A  beam  of  20-ns  laser  pulses  from  a  fixed  laser  source  was 
propagated  above  the  sea  surface  and  detected  by  a  multichannel  receiver  on  a  ship  that  traveled 
across  the  beam  line.  The  high  angular  resolution  of  the  receiver  allowed  the  direction  of  the 
source  to  be  tracked  during  the  experiment.  Representative  aerosol  distributions  were  calculated 
with  the  Advanced  Navy  Aerosol  Model  for  the  range  of  observed  weather  conditions  and  Mie 
theory  was  used  to  predict  intensity  at  the  receiver.  The  laser  beam  scattering  was  predominately 
in  the  forward  direction  and  receiver  intensity  depended  critically  on  the  beam  angle  relative  to 
the  baseline  between  source  and  receiver.  The  scattering  model  gave  reasonably  good  agreement 
with  experiment,  within  the  uncertainty  of  the  aerosol  conditions,  up  to  beam  angles  I  \fA  <  -6.5° 
where  the  scattering  was  sufficiently  strong  for  the  SBLAS  receiver.  The  receiver  also  recorded 
the  relative  arrival  times  of  the  scattered  light  for  each  of  the  angular  channels.  The  change  in 
arrival  time  with  angle,  along  with  information  about  the  source  position,  was  used  to  calculate 
the  beam  angle  if/. 
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Appendix 


The  Advanced  Navy  Aerosol  Model  (ANAM)  is  based  on  several  years  of  work  to  predict  the 
aerosol  distribution  in  a  maritime  environment  [5,7].  Here,  the  aerosols  are  considered  to  be 
uniform  spherical  particles  with  radius  r  and  refractive  index  n.  ANAM  includes  five  distinct 
types  of  aerosol  particles  each  characterized  by  a  log-normal  distribution  peaked  at  a  humidity- 
dependent  modal  radius  ps 

dN  Tp  N, 

—  =  2 L“p= - — \ — e*P 

dr  /=0V2  7z<jiexp\°i  /2 )psi 


with  a  width  parameter  o[  =  l/2m  (i  =  0,3)  and  04  =  1/2  .  Particle  0  is  a  non-hygroscopic  dust 
particle  originating  over  land  with  modal  radius  ps  =  po.  For  all  other  particles,  the  modal  radii 
increase  with  the  fractional  relative  humidity  s  (0  to  1)  according  to  ps  =  g(s)po  with  g(s)  = 
[{A -5,)/5(1^5')]1/3  Similarly,  the  complex  refractive  indices  for  particles  1-4  are  volume- 

weighted  averages  n  =  ntr+(n0-nw)g3(0)/g3(s),  where  no  is  the  index  of  the  dry  particle 
material  and  %  is  the  index  of  water  0?w  =  1.326  +  i  4.83xl0"6  at  A=  1.06  pm).  Values  of  these 
parameters  are  given  in  Table  Al.  In  this  implementation  of  ANAM,  the  total  number  density  N{ 
(cm  3)  of  each  particle  type  depends  on  the  air-mass  parameter  (amp),  wind  speed  w  (m/s)  24- 
hour  average  wind  speed  W  (m/s),  and  height  h  (m)  above  the  sea  surface, 


K-o  1 

i  -  >  amp ^  5 

[Nl=\36.55amp2] 

N0  =  0.3  x  136.55am/?2] 

„  r  amp > 5 
N,  =0.7x1 36.55am/?2  J 

N,  =  0.5462  Max [5. 866 (IT  -2.2), 0.5] 

N3  =0.007214  x  10a(0.06w)5 

N4  =  0.01x1 0A(0.07w -0.04/?). 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 
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FIGURE  CAPTIONS 


Figure  1 .  In  this  off-axis  scattering  geometry  the  laser  and  receiver  are  separated  by  R  along  the 
baseline  between  them.  The  laser  beam  is  at  an  angle  y/to  the  baseline.  Scattered  light  is  detected 
by  the  receiver  at  an  angle  6  corresponding  to  a  distance  z  along  the  beam  line  from  the  source 
and  travels  a  distance  r  to  the  receiver.  0s  is  the  angle  to  a  reference  direction  on  the  receiver. 

Figure  2.  The  laser  source  is  located  just  on  shore  at  Cowley  Beach,  Queensland,  Australia  and 
the  path  of  the  receiver  was  nearly  parallel  to  shore  and  perpendicular  to  the  beam  (red).  The 
horizontal  receiver  FOV,  shown  schematically  at  one  point  on  the  path,  was  centered  directly  to 
starboard. 

Figure  3.  Volume  scattering  functions  calculated  from  the  ANAM  cases  in  Table  1.  The 
scattering  angle  relative  to  the  beam  direction  is  6 + y/. 

Figure  4.  An  estimate  of  the  (negative)  source  angle  -0se  obtained  from  the  SBLAS  sensor  (blue 
circles)  is  compared  with  the  (negative)  source  angle  -0s  (black)  and  beam  angle  y/  (dashed) 
derived  from  the  GPS  tracking  data. 

Figure  5.  Normalized  scattered  intensity  at  the  receiver  for  three  selected  points  along  the  ship 
path  prior  to  crossing  the  beam  line.  The  model  predictions  are  calculated  on  a  0.5°  grid  with  A0 
=  3.0°,  for  the  four  aerosol  cases  of  Table  1.  The  receiver  orientations  relative  to  the  source  0s 
and  the  beam  angles  y/  were  taken  from  GPS  tracking  data. 

Figure  6.  Normalized  scattered  intensity  at  the  receiver  for  three  selected  points  along  the  ship 
path  after  crossing  the  beam  line.  The  model  predictions  are  calculated  as  in  Fig.  5. 

Figure  7.  Comparison  of  the  single-scatter  model  (black  line)  and  Monte  Carlo  calculations  (10s 
photons)  for  the  experimental  geometry  and  beam  angles  y/=  0.6°  (O),  2.4°  (□)  and  9.5°  (A)  and 
the  aerosol  scattering  function  from  Case  4.  The  MC  results  show  slightly  greater  scattering  than 
the  model  prediction  due  to  multiple  scattering  and  also  scattering  from  behind  the  beam  line. 

Figure  8.  The  relative  propagation  time  from  the  laser  to  the  receiver  determined  by  the  crossing 
time  of  the  highest  threshold  level  for  single  measurements  at  beam  angles  y/=  -6.0°  (•),  -1.5° 
(®),  1.5°  (♦),  and  6.0°  ( ^).  Timing  curves  calculated  from  (z+r)/c  using  the  known  geometry  and 
offset  to  give  the  best  fit  are  shown  as  dashed  lines. 

Figure  9.  Difference  of  beam  angle  y/c  (blue)  calculated  from  Eq.  3  using  the  estimated  source 
angle  0se  and  timing  data  <dtn/d0>  from  SBLAS  and  the  beam  angle  obtained  from  GPS  data. 
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Fig.  1.  In  this  off-axis  scattering  geometry  the  laser  and  receiver  are  separated  by  R  along  the 
baseline  between  them.  The  laser  beam  is  at  an  angle  i//to  the  baseline.  Scattered  light  is  detected 
by  the  receiver  at  an  angle  6  corresponding  to  a  distance  z  along  the  beam  line  from  the  source 
and  travels  a  distance  r  to  the  receiver.  6%  is  the  angle  to  a  reference  direction  on  the  receiver. 
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Fig.  2.  The  laser  source  is  located  just  on  shore  at  Cowley  Beach,  Queensland,  Australia  and  the 
path  of  the  receiver  was  nearly  parallel  to  shore  and  perpendicular  to  the  beam  (red).  The 
horizontal  receiver  FOV,  shown  schematically  at  one  point  on  the  path,  was  centered  directly  to 
starboard. 
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Fig.  3.  Volume  scattering  functions  calculated  from  the  AN  AM  cases  in  Table  1.  The  scattering 
angle  relative  to  the  beam  direction  is  6 + if/. 
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Fig.  4.  An  estimate  of  the  (negative)  source  angle  -Qse  obtained  from  the  SBLAS  sensor  (blue 
circles)  is  compared  with  the  (negative)  source  angle  -0s  (black)  and  beam  angle  y/  (dashed) 
derived  from  the  GPS  tracking  data. 
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Fig.  5.  Normalized  scattered  intensity  at  the  receiver  for  three  selected  points  along  the  ship  path 
prior  to  crossing  the  beam  line.  The  model  predictions  are  calculated  on  a  0.5°  grid  with  A 0  = 
3.0°,  for  the  four  aerosol  cases  of  Table  1.  The  receiver  orientations  relative  to  the  source  6>s  and 
the  beam  angles  i//  were  taken  from  GPS  tracking  data. 
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Fig.  6.  Normalized  scattered  intensity  at  the  receiver  for  three  selected  points  along  the  ship  path 
after  crossing  the  beam  line.  The  model  predictions  are  calculated  as  in  Fig.  5. 
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Receiver  Angle  O-0S  (deg) 

Fig.  7.  Comparison  of  the  single-scatter  model  (black  line)  and  Monte  Carlo  calculations  (10s 
photons)  for  the  experimental  geometry  and  beam  angles  y/=  0.6°  (O),  2.4°  (□)  and  9.5°  (A)  and 
the  aerosol  scattering  function  from  Case  4.  The  MC  results  show  slightly  greater  scattering  than 
the  model  prediction  due  to  multiple  scattering  and  also  scattering  from  behind  the  beam  line. 
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Fig.  8.  The  relative  propagation  time  from  the  laser  to  the  receiver  determined  by  the  crossing 
time  of  the  highest  threshold  level  for  single  measurements  at  beam  angles  y/=  -6.0°  (•),  -1.5° 
(®),  1.5°  (♦),  and  6.0°  (  -).  Timing  curves  calculated  from  (z+r)/c  using  the  known  geometry  and 
offset  to  give  the  best  fit  are  shown  as  dashed  lines. 
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Beam  Angle  i//  (deg) 

Fig.  9.  Difference  of  beam  angle  y/c  (blue)  calculated  from  Eq.  3  using  the  estimated  source  angle 
0se  and  timing  data  <dtu/d6>  from  SBLAS  and  the  beam  angle  obtained  from  GPS  data. 
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Table  1.  ANAM  environmental  parameters  used  to  generate  the  four  aerosol  distributions 
and  the  total  scattering  coefficient  b  and  absorption  coefficient  a  calculated  from  the  MIE 
code. 


Aerosol 

Case 

Wind  Speed 
w  (m/s) 

amp 

Height 
h  (m) 

b  (km'1) 

a  (km1) 

1 

2 

3 

2 

0.0474 

0.00023 

2 

8 

3 

2 

0.0828 

0.00058 

3 

8 

3 

300 

0.0316 

0.00006 

4 

8 

8 

2 

0.1563 

0.00219 

In  all  cases,  the  relative  humidity  was  90%  and  the  24-hour  average  wind  speed  was  3  m/s. 
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Table  Al.  ANAM  parameters  for  each  of  the  5  particle  types. 


Particle 

A)  (Mm) 

no  (A  =  1.06  pm) 

A 

B 

0 

0.03 

1.52  +  / 0.008 

- 

- 

1 

0.03 

1.392  +  / 0.000153 

1.17 

1.87 

2 

0.24 

1.47  +  /  0.0002 

1.83 

5.13 

3 

2.0 

1.47  +  /  0.0002 

1.97 

5.83 

4 

8.0 

1.47  +  /  0.0002 

1.97 

5.83 

22 


